Theory of Chiral Modulations and Fluctuations in Smectic-A Liquid Crystals Under 

an Electric Field 



O 

o 
o 

(N 
X> 

(N 



o 



Jonathan V. Sclingcr, 1 Jianling Xu, 2 Robin L. B. Selinger, 2 B. R. Ratna, 1 and R. Shashidhar 1 

1 Center for Bio/Molecular Science and Engineering, Naval Research Laboratory, Code 6900, 
4555 Overlook Avenue, SW, Washington, DC 20375 
2 Physics Department, Catholic University of America, Washington, DC 20064 

(January 24, 2000) 

Chiral liquid crystals often exhibit periodic modulations in the molecular director; in particular, 
thin films of the smectic-C* phase show a chiral striped texture. Here, we investigate whether similar 
chiral modulations can occur in the induced molecular tilt of the smectic-A phase under an applied 
electric field. Using both continuum elastic theory and lattice simulations, we find that the state 
of uniform induced tilt can become unstable when the system approaches the smectic- A-smectic- 
C* transition, or when a high electric field is applied. Beyond that instability point, the system 
develops chiral stripes in the tilt, which induce corresponding ripples in the smectic layers. The 
modulation persists up to an upper critical electric field and then disappears. Furthermore, even 
in the uniform state, the system shows chiral fluctuations, including both incipient chiral stripes 
and localized chiral vortices. We compare these predictions with observed chiral modulations and 
fluctuations in smectic-A liquid crystals. 

PACS numbers: 61.30.Gd, 64.70.Md 
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I. INTRODUCTION 

Molecular chirality leads to the formation of many 
types of modulated structures in liquid crystals |J . On a 
molecular length scale, the fundamental reason for these 
modulations is that chiral molecules do not pack paral- 
lel to their neighbors, but rather at a slight twist angle 
with respect to their neighbors. On a more macroscopic 
length scale, molecular chirality leads to a continuum free 
energy that favors a finite twist in the director field. This 
favored twist leads to bulk three-dimensional phases with 
a uniform twist in the molecular director, such as the 
cholesteric phase and the smectic-C* phase. It also leads 
to more complex phases with periodic arrays of defects, 
such as the twist-grain-boundary phases. 

In this paper, we consider the possibility of a new type 
of chiral modulation in liquid crystals. The smectic-A 
phase of chiral molecules is known to exhibit the electro- 
clinic effect: an applied electric field in the smectic layer 
plane induces a molecular tilt 0. This induced tilt is 
generally assumed to be uniform in both magnitude and 
direction. Here, we investigate whether the uniform clcc- 
troclinic effect can become unstable to the formation of a 
chiral modulation within the layer plane. There are three 
motivations for examining this possibility — one theoreti- 
cal and two experimental. 

a. Theory: The main theoretical motivation for ex- 
amining this possibility is that thin films of chiral liquid 
crystals in the smectic-C* phase show chiral modulations 
within the layer plane. These modulations have been ob- 
served in polarization micrographs of freely suspended 
films H , and have been explained using continuum elas- 
tic theory j| [?]]. In these modulations, the molecules 
form striped patterns of parallel defect walls separating 



regions with the favored chiral twist in the molecular di- 
rector. Within the narrow defect walls, the magnitude of 
the molecular tilt is different from the favored value in 
the smectic-C* phase. 

Given that in-plane chiral modulations occur in the 
smectic-C* phase, it is natural to ask whether analo- 
gous modulations can occur in the smectic-A phase under 
an applied electric field. Naively, there are two possi- 
ble answers to that question. First, we might say that 
the smectic-A phase under an electric field has the same 
structure as the smectic-C* phase, because they both 
have order in the molecular tilt. Indeed, under an electric 
field there is not necessarily a phase transition between 
smectic-A and C* |^||. Thus, we might argue that the 
smectic-A phase under an electric field should have in- 
plane chiral modulations of the form shown in Fig. 1, just 
as the smectic-C* phase does. On the other hand, the 
applied electric field itself breaks rotational symmetry in 
the smectic layer plane, and hence it favors a particular 
orientation of the molecular director. Any modulation in 
the director away from that favored orientation costs en- 
ergy. For that reason, we might argue that the smectic-A 
phase under an electric field should not have any chiral 
modulations. Because these two naive arguments contra- 
dict each other, we must do a more detailed calculation 
to determine whether chiral modulations can occur in the 
smectic-A phase under an electric field. 

b. Experiment 1: Apart from these theoretical con- 
siderations, striped modulations have been observed ex- 
perimentally in the smectic-A phase under an applied 
electric field [ [lO||Tl| ] ■ The experimental geometry is shown 
in Fig. 2(a). The chiral liquid crystal KN125 is placed in 
a narrow cell (2 to 25 /mi wide), and an electric field is 
applied across the width of the cell. One would expect 
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the smectic layers to have a uniform planar "bookshelf" 
alignment, with the layer normal aligned with the rub- 
bing direction on the front and back surfaces of the cell. 
However, the layers actually form a striped pattern, as 
shown in Fig 2(b). In the thicker cells, the striped pattern 
is quite complex, with two distinct modulations superim- 
posed on top of each other jl2|,[l3) . The main modulation 
has a wavelength of approximately twice the cell thick- 
ness, while the higher-order modulation has a wavelength 
of approximately 4 /im, regardless of cell thickness. Fur- 
thermore, the higher-order modulation is oriented at a 
skew angle of approximately 15° with respect to the main 
modulation, giving the whole pattern the appearance of 
a woven texture. 

The main modulation in these cells has been explained 
theoretically as a layer buckling instability. When an 
electric field is applied, the molecules tilt with respect 
to the smectic layers, and hence the layer thickness de- 
creases. Because the system cannot generate additional 
layers during the experimental time scale, the layers 
buckle to fill up space. This model of layer buckling pre- 
dicts layer profiles that are consistent with x-ray scat- 
tering from stripes in 2 fim cells, which do not have 
the higher-order modulation E3|. However, this model 
does not explain the observed higher-order modulation 
in thicker cells. So far, the only proposed explanation 
of the higher-order modulation has been a second layer 
buckling instability in surface regions near the front and 
back boundaries of the cell p3| . In this paper, we will 
consider whether chiral stripes can provide an alterna- 
tive explanation for this higher-order modulation. Such 
an explanation seems at least initially plausible because 
the 15° skew angle between the two observed modula- 
tions suggests a chiral effect. 

c. Experiment 2: A separate experimental motiva- 
tion for considering chiral modulations in the smectic-A 
phase comes from measurements of the circular dichroism 
(CD), the differential absorption of right- and left-handed 
circularly polarized light. Recent experiments have mea- 
sured the CD spectrum of KN125 in the smectic-A 
phase When the light propagates normal to the 

smectic layers, the CD signal is undetectable. However, 
when the light propagates in the smectic layer plane, in 
narrow cells as in Fig. 2(a), the CD signal is much larger, 
and it is quite sensitive to both electric field and temper- 
ature in a non-monotonic way. The measured CD signal 
indicates that the liquid crystals must have some chiral 
twist in the layer plane. This chiral twist might arise from 
a bulk phenomenon — either chiral modulations or chiral 
fluctuations due to incipient chiral modulations jl6| . Al- 
ternatively, it might arise from a surface phenomenon, 
such as the surface electroclinic effect [jl5[ . Hence, we 
would like to predict the bulk chiral modulations and 
fluctuations in the smectic-A phase as a function of elec- 
tric field and temperature, and assess whether such ef- 
fects can explain the CD results. 

Based on those three motivations, in this paper we 
propose a theory for chiral modulations and fluctuations 



in the smectic-A phase under an applied electric field. 
This theory uses the same type of free energy that has 
earlier been used to explain chiral stripes in the smectic- 
C* phase |-|], with modifications appropriate for the 
smectic-A phase. We investigate this model using both 
continuum elastic theory and lattice Monte Carlo simu- 
lations. In continuum elastic theory, we make an ansatz 
for the chiral modulation and insert this ansatz into the 
free energy functional. By minimizing the free energy, we 
determine whether the induced molecular tilt is uniform 
or whether it is modulated in a chiral striped texture. 
In the lattice Monte Carlo simulations, we allow the sys- 
tem to relax into its ground state, which may be cither 
uniform or modulated, without making any assumptions 
about the form of the chiral modulation. Both calcula- 
tions show that the uniform state can become unstable 
when the temperature decreases toward the smectic-A- 
smectic-C* transition, or when a high electric field is ap- 
plied. Beyond that instability point, the system develops 
a chiral modulation in the molecular tilt, and this tilt 
modulation induces a corresponding striped modulation 
in the shape of the smectic layers. This modulation is 
consistent with the higher-order stripes observed in the 
smectic-A phase. 

In addition to these theoretical results for the chiral 
modulation, we also consider chiral fluctuations that oc- 
cur before the onset of the chiral modulation itself. The 
continuum elastic theory predicts the magnitude of the 
incipient chiral stripes as a function of both electric field 
and temperature. The lattice Monte Carlo simulations 
show the incipient chiral stripes as well as localized chi- 
ral vortices in the tilt direction. These predicted fluc- 
tuations might be visible in future optical experiments. 
However, our predictions for these bulk fluctuations dif- 
fer from the CD results mentioned above |l5|] in some 
important details, including the predominant direction of 
the fluctuations and the dependence of the fluctuations 
on field and temperature. Hence, those CD experiments 
must be showing a chiral surface phenomenon, such as 
the surface electroclinic effect. 

The plan of this paper is as follows. In Sec. II, we pro- 
pose the free energy and use continuum elastic theory to 
predict chiral modulations in the tilt and the layer shape. 
In Sec. Ill, we work out the consequences of this theory 
for chiral fluctuations in the uniform phase. In Sec. IV, 
we present the lattice Monte Carlo simulations of chiral 
modulations and fluctuations in this model. Finally, in 
Sec. V, we discuss the results and compare them with 
experiments. 

II. CHIRAL MODULATIONS 

A. Tilt Modulation 

In this theory, we begin by considering a single smectic 
layer in the xy plane. Let c(x,y) be the tilt director, 
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i.e. the projection of the three-dimensional molecular 
director n(x, y) into the layer plane. The free energy can 
then be written as 



F= dA 
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Here, the r and u terms are the standard Ginzburg- 
Landau expansion of the free energy in powers of c. 
Near the transition from smectic-A to smectic-C, we have 
r — a{T — Tac)- The K$ and Kb terms represent the 
Frank free energy for splay and bend distortions of the 
director field, respectively. The b and A terms are both 
chiral terms. The b term represents the interaction of 
the applied electric field E with the molecular director, 
and the A term gives the favored variation in the direc- 
tor due to the chirality of the molecules. (This term is 
written as |c| 2 z ■ V x c rather than just z • V x c because 
the latter term is a total divergence, which integrates to 
a constant depending only on the boundary conditions. 
By contrast, |c| 2 z • V x c is not a total divergence because 
the factor of |c| 2 couples variations in the magnitude of 
c with variations in the orientation.) This free energy is 
identical to the free energy that has been used in studies 
of smectic-C* films except for two small changes. 

First, we now take the coefficient r to be positive, which 
is appropriate for smectic-A films without spontaneous 
tilt order. Second, we have added the electric field term, 
which gives induced tilt order. 

We can now ask what configuration of the tilt director 
c(x, y) minimizes the free energy. In particular, is the 
optimum tilt director uniform or modulated? To answer 
this question, we make an ansatz for c(x, y) that can de- 
scribe both the uniform and modulated states. Suppose 
the electric field is in the y direction, E = Ey, which 
favors a tilt in the x direction. In a uniform state, the 
system has the electroclinic tilt c(x, y) = cox. By com- 
parison, in a modulated state, the system has a chiral 
striped pattern of the form shown in Fig. 1. Here, the 
tilt director is modulated about the average value of cq 
in the x direction. The magnitude of the tilt is larger 
when the tilt varies in the direction favored by molecular 
chirality, and it is smaller when the tilt varies in the oppo- 
site direction. This chiral modulation can be represented 
mathematically by the ansatz 



c x (x,y) = c + ci cos(q x x + q y y), 
c y (x,y) = ci sin(q x x + q y y). 



(2a) 
(2b) 



This ansatz has four variational parameters: cq gives the 
average tilt, c\ gives the amplitude of the modulation, 
and q x and q v give the wavevector. Within this ansatz, 
c\ = corresponds to the uniform state, and c± ^ cor- 
responds to the modulated state. This is not the most 
general possible ansatz. In general, the c x and c y compo- 
nents of the modulation might have different amplitudes, 



and might not be exactly 90° out of phase. Furthermore, 
the modulation might have multiple Fourier components. 
These possibilities will be considered in the lattice Monte 
Carlo simulations of Sec. IV. For now, this simple ansatz 
demonstrates the relevant physics of the modulation. 

To determine whether the optimum tilt director is uni- 
form or modulated, we insert the ansatz into the free en- 
ergy and average over position. The resulting free energy 
per unit area can be written as 
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(3) 



Here, K — ^(Ks + Kb) is the mean Frank constant. 
Minimizing the free energy over q x and q y gives 



q x = 



Acq 
K ' 



q y = 0. 



(4a) 
(4b) 



These expressions give the wavevector of the first unsta- 
ble mode of chiral modulation. Note that this wavevector 
is in the x direction. This result is reasonable, because 
the chiral stripes in the smectic-C* phase also have the 
modulation wavevector parallel to the average tilt direc- 
tion. Inserting these expressions back into the free energy 
gives the result 
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(5) 



expressed in terms of the two parameters cq and c\. 

Note that this free energy is only thermodynamically 
stable (bounded from below) for a certain range of A. 
The combination of the u and A terms can be regarded 
as a quadratic form in the variables Cq and c 2 . The free 
energy is thermodynamically stable only if this quadratic 
form is positive-definite, which requires A 2 < 3Ku. If 
this condition is not satisfied, then the free energy must 
be stabilized either by higher-order terms (such as |c| 6 ) 
not considered here, or by the constraint |c| < 1, which 
results from the fact that c is the projection of the three- 
dimensional molecular director n into the layer plane. 

We can now find the minimum of the free energy with 
respect to cq and c\. The extrema of the free energy are 
given by 

d F ( A 2 A 

a = rc ° + uc ° + 2 ( u ~ 2kj C ° C ^ ~ bE = °' ^ 



d F 
dc^A 
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2K 



CqCi = 0. 



(6b) 



One solution of these equations is the uniform electro- 
clinic state, in which c\ — and Co is given implicitly 
by 
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rc + ucl - bE = 0. 



(7) 



The question now is whether this solution is a minimum 
or a saddle point of the free energy. If it is a minimum, 
the uniform state is stable; if it is a saddle point, the uni- 
form state is unstable to the formation of a chiral modu- 
lation. To answer this question, we calculate the Hessian 
matrix of second derivatives at the uniform solution. The 
determinant of this matrix is 
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(8) 



If the chiral coefficient A is weak, with A 2 < 2Ku, then 
the determinant is always positive, implying that the 
uniform state is always stable. On the other hand, if 
A 2 > 2Ku, then the determinant can become negative, 
and the uniform state can become unstable. In that case, 
the system has a second-order transition from the uni- 
form state to the modulated state when the determinant 
passes through zero. This transition can be driven by 
increasing the electric field, which increases cq. It can 
also be driven by decreasing the temperature toward the 
smectic-A-smectic-C* transition, which reduces r and in- 
creases Co- 

By combining Eqs. (^-^), we can calculate properties of 
the uniform and modulated states around the transition. 
First, the transition occurs at the electric field E* given 

by 



a-.i-W-^l* ^f^-il. (9) 
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On the uniform side of the transition, the tilt is 
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(10) 



On the modulated side of the transition, the amplitude 
of the chiral modulation increases as 
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Furthermore, the wavevector of the chiral modulation is 



AcS 



A /r\i/2 / A 



K K \u 



Ku 



-1/2 



(12) 



on the modulated side of the transition. 

There is one more constraint on the uniform- 
modulated transition that we have not considered yet. 
Because the tilt director c is the projection of the three- 
dimensional molecular director n into the layer plane, its 
magnitude is limited to |c| < 1. This constraint implies 



that there is an upper critical field E msix beyond which 
the director is locked at unit magnitude in the x direction. 
From Eq. (Q) , the upper critical field can be estimated as 



En 



b 1 (r + u). 



(13) 



Above this field, the chiral modulation disappears and 
the system becomes uniform again. 

By combining the results derived above, we obtain a 
phase diagram in terms of the electric field E and the 
chiral coefficient A, which is shown in Fig. 3. For large 
values of A, the system has a transition from the uniform 
state to the modulated state at the field E* , and then it 
goes back into the uniform state at the field E max . At 
A 2 w K(r + 2u), these two transitions intersect. Below 
that value of A, the system remains in the uniform state 
for all electric fields, with no modulated state. 



B. Layer Modulation 

So far, we have considered only modulations in the tilt 
director in a flat smectic layer. However, a recent theory 
of the Pp' rippled phases of lipid membranes shows that 
any modulation in the tilt director will induce a mod- 
ulation in the curvature of the membrane jl7|,[l8). The 
same theoretical considerations that apply to modula- 
tions in lipid membranes also apply to smectic layers in 
thermotropic liquid crystals JTs| ] - Hence, we can carry 
over these theoretical results to predict the curvature 
modulation induced by the chiral stripes in the smectic- 
A phase. Let h be the height of the smectic layer above 
a flat reference plane. The results of Refs. |l7],[l^] then 
predict 



cPh 
dx 2 



7 dc x 
k dx 



A 



HP 



(14) 



for a modulation in the x direction. Here, k is the cur- 
vature modulus of the layer, 7 is the non-chiral coupling 
between tilt and curvature, and Ahp is the chiral cou- 
pling [^0). After inserting the tilt modulation of Eq. (||), 
with q y = as found above, we can integrate this differ- 
ential equation to obtain the curvature modulation 



h(x) 



7£i _ AhpCqCi 

nq x nql 



suvq x x 



^HPC? . 

— — T sin 2q x x. 



(15) 



The most important feature to notice about this modu- 
lation is that it includes both sinq^a; and sin 2q x x terms. 
As a result, the layer profile has the shape shown in 
Fig. 4(a) (with a highly exaggerated amplitude). This 
modulation has the symmetry — it has a two- fold 
rotational symmetry, but it does not have a reflection 
symmetry in the xy or yz plane because of the chiral 
coupling Ahp- 
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Although our model applies only to a single smectic 
layer, the symmetry of the modulation gives some in- 
formation about the packing of multiple smectic layers. 
Because a single layer does not have a reflection symme- 
try in the xy or yz plane, the packing of multiple lay- 
ers should not have such a symmetry either. Instead, 
the most efficient packing of multiple layers should have 
the form shown in Fig. 4(b), with a series of sawtooth- 
type waves stacked obliquely on top of each other. For 
that reason, this chiral instability should lead to oblique 
stripes, which are not parallel to the average layer nor- 
mal along the z axis. The three-dimensional wavevector 
therefore has a q z component as well as a q x compo- 
nent. This inclination of the stripes gives a macroscopic 
manifestation of the chiral mechanism that generates the 
stripes. 



III. CHIRAL FLUCTUATIONS 

At this point, let us return to the problem of tilt varia- 
tions in a single smectic layer. In the previous section, we 
showed that a chiral instability can give a periodic mod- 
ulation in the molecular tilt. Even if the system does not 
have a periodic chiral modulation, it can still have chiral 
fluctuations about a uniform ground state, i.e. ci-type 
fluctuations about the uniform electroclinic tilt cq. Such 
fluctuations can have an important effect on the optical 
properties of the system. For that reason, in this sec- 
tion we investigate the theoretical predictions for chiral 
fluctuations in the uniform phase. 

The first issue in this calculation is to determine what 
quantity gives an appropriate measure of the strength of 
chiral fluctuations. The simplest possibility is just the 
chiral contribution to the free energy of Eq. (|J) , 



chiral 



-XqxCo^. 



(16) 



An alternative possibility is suggested by a recent model 
for the transition between the isotropic phase and the 
blue phase III pi] ]. This work proposed the chiral order 
parameter 



tp = Q ■ V x Q = eijkQudjQ 



kl; 



(17) 



where Qij = riiUj — \b~ij is the tensor representing ne- 
matic order. In our problem, the three-dimensional ne- 
matic director n can be written as 



(l-||cf)z, 



and hence the tensor becomes 



2 1 

c x C X Cy c x 

2 1 

C X Cy C y — g Cy 

C X Cy g — C x — C y 



(18) 



(19) 



Into this expression we insert the ansatz of Eq. (g) for the 
fluctuations, with q y = 0. After averaging over position, 
the chiral order parameter simplifies to 



(20) 



Remarkably, this result for the chiral order parameter is 
equivalent to the expression for Fdiirai, up to a constant 
factor. This equivalence shows that either expression can 
be used as a measure of the strength of chiral fluctuations. 

We can now calculate the expectation value of the chi- 
ral fluctuations. Applying the equipartition theorem to 
the free energy of Eq. @ gives 



k B T\ 2 c 2 



chiral 



Kr - (A 2 - 2uK)cl 



(21) 



For low fields, we have cq = bE/r. Furthermore, we have 
r = a(T — Tac) near the smectic- A-smectic-C transi- 
tion. Under those circumstances, the expectation value 
becomes 



k B T\ 2 b 2 E 2 



al 



Ka 3 (T ~ Tac) 3 - (A 2 - 2uK)b 2 E 2 ' 



(22) 



This expression shows that the system has chiral fluctu- 
ations in the uniform state. The magnitude of the chiral 
fluctuations depends on both the applied electric field 
and the temperature. In particular, this measure of the 
chiral fluctuations begins at (|F c hirai|) = for zero field 
and increases as the applied electric field increases. If 
A 2 > 2Ku, then increasing the field drives the system to- 
ward the uniform-modulated transition, where the mag- 
nitude of the fluctuations diverges. Otherwise, increas- 
ing the field drives the system toward a finite asymptotic 
value of the chiral fluctuations. 



IV. MONTE CARLO SIMULATIONS 

To investigate further this model for chiral modula- 
tions and fluctuations, we have done a series of Monte 
Carlo simulations. These simulations serve two main 
purposes. First, they allow the system to relax into its 
ground state, which may be either uniform or modulated, 
without the need for any assumptions about the form of 
the chiral modulation. Hence, they provide a test of the 
assumed form of the chiral modulation in Eq. (||). Sec- 
ond, the simulations provide snapshots of the tilt director 
field for different values of the electric field E and the chi- 
ral coefficient A. Thus, they help in the visualization of 
the chiral modulations and fluctuations. 

In the simulations, we represent the tilt director in a 
single smectic layer by a discretized xy model on a two- 
dimensional hexagonal lattice. Each lattice site i has 
a tilt director Ci of variable magnitude |ci| < 1. We 
suppose that K$ = Kb = K. The discretized version of 
the free energy of Eq. ([!]) then becomes 
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(23) 



Here, is the unit vector between neighboring sites i 
and j. The factor of | is required because of the hexag- 
onal coordination of the lattice. This model is similar to 
a discretized model for tilt modulations studied in the 
context of Langmuir monolayers p2| . 

We use a lattice of 100 x 100 sites with periodic bound- 
ary conditions. We fix the parameters r — 0.5, u = 1, 
6=1, and K = 1.5, and vary the chiral coefficient A and 
the strength of the electric field E in the y direction. For 
each set of A and E, we begin the simulations at a high 
temperature, with all the directors Cj = 0, corresponding 
to an untilted smectic-A phase. We allow the system to 
come to equilibrium and then slowly reduce the temper- 
ature, so that the system can find its ground state. This 
procedure can be regarded as a simulated-annealing min- 
imization of the discretized free energy with a fixed set of 
parameters. To refine the phase diagram further, some 
simulations of the modulated phase were performed using 
system size 100 x 4. The results are in close agreement 
with those in the 100 x 100 system, since the modulation 
is essentially one-dimensional. 

Figure 5 shows the director configuration that forms 
in a typical run with A = 2.4 and E = 0.4. Here, the sys- 
tem has relaxed into a configuration of chiral stripes with 
wavevector in the x direction, perpendicular to the elec- 
tric field applied in the y direction. This configuration 
appears similar to the ansatz proposed in Fig. 1. The 
simulation results can be compared quantitatively with 
the predictions of continuum elastic theory. For these 
parameters, Eq. (||) predicts E* « 0.4 and Eq. ( p"3| ) pre- 
dicts -Emax = 1.5. Hence, the system is quite close to the 
uniform-modulated transition at E = E* and far from 
-Emax- From Eqs. (|o|) and (|l2|), the average value of the 
tilt should be Cg ss 0.5 and the modulation wavevector 
should be q* » 0.8, corresponding to a wavelength of 
27r/<7* w 8 lattice units. These predictions are approxi- 
mately consistent with the simulation results. 

For a further comparison of the continuum elastic the- 
ory with the simulations, we can look at the shape of 
Cy(x) as a function of x. The ansatz of Eq. (^|) assumes 
that this is a sine wave with a single wavevector q x . In 
the simulation, the modulation can take any shape, so we 
can assess whether the shape is actually sinusoidal. Fig- 
ure 6 shows the shape of the modulation for two sets of 
parameters. In Fig. 6(a) we have A = 2.25 and E — 1.48, 
which is quite close to £ max w 1.5. For these parameters, 
the amplitude of the modulation is very small, and the 
shape of the modulation is very well fit by a sine wave. 
By contrast, in Fig. 6(b) we have A = 2.55 and E = 0.2, 
which is far from the transitions at E* and E max . Here, 
the amplitude of the modulation is much larger, and it is 
clearly not sinusoidal. This comparison shows that the 



theoretical assumption of a sinusoidal modulation is ap- 
propriate close to the transitions, where the modulation 
amplitude is small, but it is not appropriate well inside 
the modulated state, where the amplitude is large. 

The simulation results for the phase diagram are sum- 
marized in Fig. 7. Note that this phase diagram has the 
same structure as the phase diagram predicted by con- 
tinuum elastic theory in Fig. 3. For large A, the system 
goes from the uniform state to the modulated state at 
the electric field E* . The system remains in the modu- 
lated state up to the field -Emax, at which point it goes 
back into the uniform state. As A decreases, the range 
of the modulated state in the phase diagram decreases, 
and it finally vanishes at A « 2.14. For smaller values 
of A, the system stays in the uniform state for all values 
of the electric field. The numerical value of the phase 
boundary E*(X) is shifted somewhat from the theoreti- 
cal prediction. The difference can be attributed to the 
discretization of the system, which changes the free en- 
ergy of the modulated state by a few percent. This small 
change in the free energy is enough to give a noticeable 
shift in E*(X). 

In addition to these results for the phase diagram in 
the thermodynamic ground state, the Monte Carlo simu- 
lations also allow us to observe chiral fluctuations in this 
model. As noted above, the simulation procedure in- 
volves beginning in a high-tcmperature disordered state, 
allowing the system to come to equilibrium, and then 
gradually reducing the temperature to zero. If we inter- 
rupt the simulations at a nonzero temperature, before 
the system reaches the ground state, then we can take a 
snapshot of the fluctuations. For example, Fig. 8 shows 
a snapshot of the simulations at A = 2.4 and E = 0.2, 
which has been interrupted at temperature T — 0.081. 
In the thermodynamic ground state given by the phase 
diagram, this system is uniform. At this finite tempera- 
ture, the uniform state shows chiral fluctuations, which 
take the form of incipient chiral stripes. These incipient 
chiral stripes show a specific realization of the chiral fluc- 
tuations discussed in Sec. III. As the field is increased, 
these fluctuations grow in magnitude and eventually be- 
come equilibrium stripes at E*(\). 

In addition to the incipient chiral stripes, the model 
also shows another type of chiral fluctuations, which are 
localized chiral vortices in the tilt director. Figure 9 
shows an example of the vortices for A = 2.4 and E = 0. 
These vortices seem to be forming an incipient hexagonal 
lattice. The vortex pattern is a nonequilibrium fluctua- 
tion, which anneals into the uniform state as the simu- 
lation temperature is decreased. When a small electric 
field is applied, the vortices are suppressed. The field 
strength required to suppress vortices increases as A in- 
creases. This nonequilibrium vortex pattern is similar 
to the equilibrium hexagonal lattice of vortices that has 
been predicted in theoretical studies of the chiral smectic- 
C* phase The nonequilibrium vortex pattern prob- 

ably evolves into the equilibrium vortex lattice as r is 
decreased from the smectic-A phase into the smectic-C* 
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phase, but we have not yet tested this scenario in the 
simulations. 



V. DISCUSSION 

In the preceding sections, we have answered the the- 
oretical question that motivated this study. Our model 
shows that the smectic-A phase under an applied elec- 
tric field can become unstable to the formation of a chi- 
ral modulation within the layer plane, which is similar 
to the chiral striped modulation that has been observed 
in thin films of the smectic-C* phase. The transition 
from the uniform to the modulated state occurs when 
the temperature decreases toward Tac or when a high 
electric field is applied. It is somewhat surprising that 
an applied electric field can induce a modulation in the 
director away from the orientation favored by the field. 
However, an analogous effect has been predicted by a 
recent study of cholesteric liquid crystals in a field p3[ . 
That theoretical study showed that a high electric field 
can induce a transition from a paranematic phase to a 
cholesteric phase, in which the director has a modulation 
away from the field direction. In both that problem and 
our current problem, the transition between the uniform 
and modulated states is controlled by the competition 
between field-induced alignment and chiral variations in 
the orientation of the field- induced order parameter. 

We can now compare this model with the two experi- 
mental results mentioned at the beginning of this paper. 
In the first experiment, two types of stripes were observed 
in the smectic-A phase of the chiral liquid crystal KN125 
under an applied electric field — a main modulation with 
a wavelength of approximately twice the cell thickness 
and a higher-order modulation with a wavelength of ap- 
proximately 4 /im, independent of cell thickness |^2|^3[ . 
The main modulation has been explained as a layer buck- 
ling instability. Our model of chiral stripes provides 
a possible explanation of the higher-order modulation. 
In particular, it predicts the observed symmetry of the 
modulation — the observed skew angle between the main 
stripes and the higher-order stripes in Fig. 2(b) corre- 
sponds to the theoretical packing angle in Fig. 4(b). Fur- 
thermore, it predicts that the wavelength of the higher- 
order modulation is determined by material properties of 
the liquid crystal, not just by the cell thickness. 

One possible objection to this model for the exper- 
iment is that the observed stripe wavelength does not 
depend on electric field, while the prediction of Eq. (^a) 
depends on electric field implicitly through the uniform 
tilt cq. The response to this objection is that the stripe 
wavelength is not necessarily in equilibrium. Rather, the 
stripe wavelength is probably determined by the wave- 
length at the onset of the instability and cannot change 
in response to further variations in electric field. A sec- 
ond possible objection is that the predicted stripes occur 
only close to Tac, while the observed stripes occur in 



the liquid crystal KN125, which does not have a smectic- 
A-smectic-C* transition. The response to that objec- 
tion is that KN125 is an unusual liquid crystal with a 
large electroclinic effect over a surprisingly wide range 
of temperature p3j. For that reason, this material can 
show chiral stripes over a wide range of temperature. A 
further test of the theory would be to see whether the 
higher-order stripes occur in a liquid crystal that has the 
standard temperature-dependent electroclinic effect near 
Tac, an d to see whether these stripes are more sensitive 
to temperature. 

The second experiment mentioned at the beginning of 
this paper measured the CD spectrum of KN125 in the 
smectic-A phase under an applied electric field [Tq j. A 
very large CD signal was found for light propagating in 
the smectic layer plane, much larger than the CD signal 
for light propagating normal to the smectic layers. This 
large CD signal indicates that the liquid crystals have 
some chiral twist in the smectic layer plane. This twist 
might arise from a bulk phenomenon, such as the chiral 
fluctuations of the uniform smectic-A phase considered 
in this paper. Alternatively, it might arise from a surface 
phenomenon, such as the surface electroclinic effect. 

Although our model for chiral fluctuations in the uni- 
form smectic-A phase gives one possible source for a CD 
signal that could be measured in optical experiments, 
our predictions differ from the experimental results in 
two important details. First, the predominant wavevec- 
tor of the predicted chiral fluctuations is in the x di- 
rection, perpendicular to the electric field. By contrast, 
the experiments are sensitive to chiral fluctuations in the 
y direction, along the electric field, because the light is 
propagating in that direction. Second, the theory pre- 
dicts that the quantity (|-Fchiral|)> measuring the strength 
of chiral fluctuations, should be zero at zero field and 
should increase monotonically with increasing field. In 
the experiments, the measured CD signal is nonzero at 
zero field, and it can increase or decrease with increasing 
field, depending on the temperature. Thus, the experi- 
ment must be showing a chiral surface phenomenon. A 
specific model based on the surface electroclinic effect has 
been presented in Ref. [|l5| . The bulk chiral fluctuations 
discussed in this paper might be observable in future op- 
tical experiments, especially if the surface effects can be 
suppressed through appropriate surface treatments. 

In conclusion, we have shown that the uniform electro- 
clinic effect in the smectic-A phase of chiral liquid crystals 
can become unstable to the formation of a chiral modu- 
lation in the layer plane. In the modulated state, there 
are stripes in the orientation of the molecular director, 
analogous to the stripes that have been observed exper- 
imentally in thin films of the smectic-C* phase. The 
same mechanism also gives chiral fluctuations in the uni- 
form smectic-A phase, which grow in magnitude as the 
modulated state is approached. These theoretical results 
provide a possible explanation for stripes observed in the 
smectic-A phase, and they predict chiral fluctuations that 
may be observed in future optical experiments. 
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FIG. 1. Ansatz for the chiral modulation of the molecular 
tilt in a layer of the smectic-A phase under an applied electric 
field. The arrows represent the projection of the molecular 
tilt into the layer plane. The electric field is applied in the y 
direction, and the average molecular tilt is in the x direction. 
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FIG. 2. (a) Experimental geometry of the striped modula- 
tion in the smectic-A phase in narrow cells (from Ref. p^|). 
(b) Polarization micrograph of the striped pattern in a 15 /im 
cell, showing the xz plane (with the line of sight along the 
y axis). Note that there are two distinct modulations with 
wavelengths of approximately 30 /im and 4 fim, oriented at an 
angle of approximately 15° with respect to each other (from 
Ref. @). 
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FIG. 3. Theoretical phase diagram in terms of A and E 
for fixed r = 0.5, u = 1, 6 = 1, and K = 1.5, showing the 
uniform- modulated transitions at E*(\) and E ulax (X). 
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FIG. 4. (a) Profile of the height modulation of a single 
smectic layer, showing the Cj symmetry. The vertical dis- 
placement is highly exaggerated, (b) The most efficient pack- 
ing of multiple smectic layers with a modulation of this sym- 
metry. The dashed lines indicate the orientation of the stripes 
in the xz plane. 



FIG. 5. Simulation results for the molecular tilt in a layer 
of the smectic-A phase, showing a chiral modulation for the 
parameters r = 0.5, u = 1, b = 1, K = 1.5, A = 2.4, and 
E — 0.4. The lines represent the projection of the molecular 
tilt into the layer plane. The applied electric field is in the y 
direction. 




FIG. 6. Simulation results for the shape of c y (x): (a) Near 
a transition, for A = 2.25 and E = 1.48. (b) Well inside the 
modulated state, for A = 2.55 and E — 0.2. Note that the 
modulation has a sinusoidal form in the first case, but not in 
the second. 
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X 

FIG. 7. Simulation results for the phase diagram in terms 
of A and E for fixed r = 0.5, u = 1, b = 1, and A' = 1.5. FIG. 9. Finite-temperature simulation results for the pa- 

The symbols indicate the observed transitions between the rameters r = 0.5, u = 1, b = 1, K = 1.5, A = 2.4, and E = 0, 
uniform and modulated states. The discrepancy between this showing localized chiral vortices, 
figure and Fig. 3 is due to the discretization of the lattice. 
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FIG. 8. Finite-temperature simulation results for the pa- 
rameters r = 0.5, u = 1, b = 1, K = 1.5, A = 2.4, and 
B = 0.2, showing incipient chiral stripes. 
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